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Abstract 

We analyze how the thermal history of the universe is influenced by the statistical description, 
assuming a deviation from the usual Bose-Einstein, Fermi-Dirac and Boltzmann-Gibbs distribution 
functions. These deviations represent the possible appearance of non-extensive effects related with the 
existence of long range forces, memory effects, or evolution in fractal or multi-fractal space. In the early 
universe, it is usually assumed that the distribution functions are the standard ones. Then, considering 
the evolution in a larger theoretical framework will allow to test this assumption and to place limits 
to the range of its validity. The corrections obtained will change with temperature, and consequently, 
the bounds on the possible amount of non-extensivity will also change with time. We generalize results 
which can be used in other contexts as well, as the Boltzmann equation and the Saha law, and provide 
an estimate on how known cosmological bounds on the masses of neutrinos are modified by a change in 
the statistics. We particularly analyze here the recombination epoch, making explicit use of the chemical 
potentials involved in order to attain the necessary corrections. All these results constitute the basic 
tools needed for placing bounds on the amount of non-extensivity that could be present at different eras 
and will be later used to study primordial nucleosynthesis. 
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1 Introduction 



This paper is the first in a series that will thoroughly analyze the influence of the statistical description 
in the standard picture we have for the evolution of the early universe. To do so, we shall insert the 
usual statistical mechanics, i.e. Boltzmann-Gibbs', in a larger framework, given by non-extensive theories. 
These theories are parameterized by a new degree of freedom, related to the amount of non-extensivity 
present in the system under consideration. The standard statistical mechanics is the only one that respects 
the property of extensivity, i.e. the entropy of a system formed by several subsystems equals the sum of 
the entropies of each separate subsystem. Then, considering this case as a particular situation out of a 
larger set of possibilities, our aim will be twofold. On the one hand, we shall comment on the reasons by 
which the universe as a whole could deviate from being an extensive system, and on the possible amount 
of these deviations. On the other, we shall systematically analyze how many of the standard results of the 
usual cosmological model, the hot big bang theory, get affected by such an slight change. This change will 
manifest itself in the form of the quantum distribution functions, that we shall modify at the beginning 
and use thorough the paper. Using these generalized predictions, we shall be in position to assess how 
our present knowledge of observational data can bound the statistics operative at different eras of cosmic 
evolution. The corrections obtained will change with temperature, and consequently, the bounds on the 
possible amount of non-extensivity will also change with time. 

We shall of course draw upon some previous results, for instance, the modification that a change in 
the statistics introduce in the energy density of the universe was previously studied by several authors. 
Particularly, some of us have been working in the topic of nucleosynthesis, under what it is known as the 
asymptotic approach of quantum distribution functions [[]], Q. However, we shall re-derive some of these 
results in order to have all of them within the same non-extensive framework; this will allow to consistently 
use our new results in the follow up applications. Indeed, these preliminary works will be used as the 
launch pad for a consistent analysis on the influence of statistical mechanics on all standard -textbook- 
cosmology: We shall generalize and obtain corrections due to this new setting (up to first order in the 
deviation parameter) for processes going from the decoupling of hot and cold relic species, nucleosynthesis, 
recombination, to matter-radiation equality. In our way up, we shall as well derive results which can be 
used in other contexts, as the generalized Boltzmann equation and the Saha law. A direct application of 
them is also carried on here, particularly in the section concerning neutral hydrogen formation. 

Many are the papers on astrophysics and cosmology using non-standard statistical mechanics. Among 
them, we would like to mention here the works related with the solar core ||, where an interesting analysis 
of neutrino production in a non-extensive setting was presented. It is both interesting and instructive to 
follow the recent discussion on the need for a new statistical description in the solar core, to that end see 
the papers of Refs. 0, p|. Another astrophysical application has been on the topic of high and ultra-high 
energy cosmic rays ||. In particular, it was recently proven that a non-extensive setting can not solve the 
GZK cutoff problem |7j|. Of course we must recall other previous papers on cosmology ]|], and make a 
special remark on the proof of the T oc relationship |J, the analysis of COBE satellite data as a bound 
for non-extensivity [ 10 1 , the use of precision cosmology (MAP, Planck, SDSS) with the same aim [11|, and 
some previous ideas on kinetic theory This series of papers encompass most of these results, enlarge 
them, and tries to form an unified picture of non-extensive cosmology, both bounding its possible range of 
validity and extent and, from a formal point of view, studying how a generally non-stated assumption can 
hide lots of implications. 

The layout of this paper is as follows. In the sake of making a self-contained work, especially for 
cosmologists, we shall state in Section 2 some of the most important ideas and the basic development of 
non-extensive statistics. For the ease of the discussion, some brief historical comments will be made also 
there. We shall also present the new quantum distribution functions, that we shall use in the rest of the 
work, following in this regard the discussion made by Buyiikkihg et al. and Tirnakh et al. [13]. Section 
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2, then, does not pretend to show new results. Section 3 discusses the values of the plasma parameter 
along cosmic history, and assess the a priori possibility that the universe as a whole can be considered 
as a non-extensive system. Section 4 and 5 are the basis for all further development; they state the 
thermodynamical picture and the cosmologically conserved quantities in this new setting. Then, some 
new applications follow: Section 6 is concerned with the decoupling of relativistic and non-relativistic 
species; Section 7, with the process of recombination; and Section 8, with the new Boltzmann equation 
and the study of the process of freezing. Then, a brief Section 9 states the predictions for the current 
values of some observables of cosmological importance; and Section 10 studies hot and cold relics in the 
new framework and generalizes some of the well known bounds on neutrino masses and on other species. 
Finally, in Section 11, it is briefly accounted for the correction that non-extensivity would introduce to the 
value of matter-radiation time and temperature. We end by summarizing and obtaining some conclusions 
of a general nature. Paper II in this series will be an analytical analysis of the nucleosynthesis process, 
going all the way up to the formation of Deuterium and Helium 4. We shall be specially concerned with 
the analysis of the nuclear principle of detailed balance and the influence it has in this new setting. 



2 Basics of non-extensive statistics 

Among physicist and astronomers, there is a non-stated consensus on that the Boltzmann-Gibbs (BG) 
statistical mechanics is always applicable. However, in analogy with Newtonian dynamics and Special 
Relativity, it could a priori justified to consider it just as part of a bigger framework, an enlarged statistical 
description where extensive as well as non-extensive phenomena could be taken into account fi"4}| . An 
interesting generalization of the BG entropy form has been recently proposed by Tsallis |l5| (for recent 
[1^, 17], for a full bibliography see [ 18 1 ) . This new entropy, that we discuss below, possesses the 
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usual properties of positivity, concavity and irreversibility, and generalizes the additivity in a non-extensive 
way. Examples in which BG statistics seems to present serious problems are systems for which there are 
long-range forces, and/or which present memory effects, and/or are subject to evolution in a non-Euclidean 
space [19]. 



The generalized form for the entropy is 

5 = k ^f^, q G R, (1) 

q-1 

where k a positive constant, W the total number of physical states accessible to the system, and the set of 
probabilities pi satisfies 

w 

5> = 1 - ( 2 ) 

Eq. (|l|) recovers the usual (BG) form for the entropy in the limit q — > 1, i.e. 

w 

lim S q = -k^pilnpi. (3) 

9-> i=i 

The entropic index q, intimately related to, and determined by, microscopic dynamics, characterizes the 
system under consideration. This reflects itself in the pseudo-additivity law for the entropy S q , 

S q (A + B)/k = S q {A)/k + S q {B)/k + (1 - q)[S q (A)/k] [S q (B)/k], (4) 

where A y B are two independent systems in the sense that the joint probabilities (those corresponding to 
the system A + B) are such that Pij(A + B) = pi(A)pj(B). 
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2.1 Expectation values 

We now introduce the following non-normalized expectation value, 

w 

(A) q = J2pl A i (5) 

such that (A)i corresponds to the standard mean value for the observable A. If the system is a quantum 
one, its description is given in terms of the density operator, p, with eigenvalues {pi}. Then, the generalized 
entropy is 

1 _ Trp q 

Sq = k——f- (Trp=l), (6) 

and the expectation value becomes 

(A) q = Trp q A. (7) 

Eq. can be recast as 

S q = -k(ln q p) q . (8) 

If the system is classic, and the relevant variables are continuum ones, we can describe it by a probabil- 
ity distribution p(r), where r is a dimensionless variable, say, in a many-body phase space. Then, the 
generalized entropy will be 

= k l-jdr[p{r)f w . th r ^ = ^ 



q-l 

and the expectation value, 

{A), = I df[p(r)] q A(f). (10) 



We shall present the formalism in the case in which the system is described by a set of microscopic 
probabilities, W. 

2.2 Canonical Ensemble 

The first non-trivial physical situation is that of a system in thermal contact with a thermostat at a 
temperature T. We must then extremize the entropy taking into account several constraints. The first one 
is just the definition of probability, 

w 

Ek = 1- (11) 
1=1 

For the rest, a detailed discussion is in order. 
2.2.1 Internal energy 

There is some freedom concerning the choice of the constraint imposing the relationship among the different 



energy levels with the total internal energy. The first choice was introduced by Tsallis [15]: to Eq. (11), it 
is added the following constraint 

W 



E^ = [/(1) - (i 2 ) 

i=i 
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Here, an index (1) refers to the first choice, and {e{\ are the Hamiltonian eigenvalues of the system. Using 
standard techniques, it can be seen that the {pj}-values that extremize S q with the imposed constraints 
are 

It must be said that [3* is not the Lagrange multiplier associated with the constraint upon the internal 
energy. This expression recovers the usual one (pi cx e~^ e ') in the limit q — > 1 and it depends on the 
microscopic energies as a power law instead of the familiar exponential function. 
The second choice |l5[ postulates that 

W 



E^ = c/(2) - (14) 



=i 



The {pi}-values that now extremize S q are 



pf ] = l - "~ rgP , (15) 

where we have defined the generalized partition function 

w 

4 2) = E[i-(i-^il 1/M - (16) 

This result differs from the previous in that the role played by (1 — q) is now equivalent to what was played 
by (q — 1), and that in this case, (3 is the Lagrange multiplier associated with the constraint on the internal 
energy. The probability distribution can be conveniently recast as 

pf ] = e ^Z with (1 7 ) 

Zq j=l 

what allows to prove the following series of equalities |2C], 

T dUf> 

- U?)-TS q = -±\n q Zf), (19) 
^ = "^f^, (20) 

C ( 2) = T ds q _du^ __ T ^f 

with F q , U q and C q standing for the corresponding generalizations of Helmholtz free energy, internal energy 
and specific heat. Then, the formal thermodynamical structure (Legendre transformations) remains valid. 
Three unwanted consequences, however, should be noted 21] 



[1 - (1 - (O/feip/a-*) 



1. The distribution given by Eqs. (15) and ( jig ) is not invariant against an uniform shift of the energy 
zero. 

2. The mean value of a constant differs from the constant itself, {l) q ^1. 
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3. Finally, if two systems A and B are such that pfj +B = pfpf and e^f~ a = ef + ej*, then 

Up (A + B)/k = uM(A)/k + U®{B)/k + (1 - g)[t^ 2 )(A)5 fl (5)/fc][£^)(S)5,(il)/fc] f (22) 
what differs from U q (A)/k + U q (B)/k. This means that the energy is not an additive quantity. 



A+B 



A third choice for the energy constraint was proposed by Plastino, Mendes and Tsallis [21|, and its key 
aspect is normalization. The constraint is now, 



V 14 " r> q 



(23) 



i.e. it weights the eigenvalues with a set of escort probabilities p\IY^i=\P\- The set {pi} which now 
extremizes S q are (see the Appendix for the definition of the function involved) 



P. 



(3) 



[1 - (1 - q)(3(ei - U^)/EZi (P? } ) q ] m - q) _ ex P J-^(6, - U^/ZZi (p?) q ] 



Z 



(3) 



*(3) 



(24) 



where 



W 

= e 

3=1 



l-(l-«))9(c i -Efi)/E(W 3 T 



W n V(l-9) 

(3)n 



t=l 



VP 

E eX P« 

i=i 



-/%-^ (3) )/£(pi 



(3). 



8=1 



(25) 



It can be shown that, if T = l/k/3 [21], 



1 

T 



77(3) 
1 <1 



dS n 



(3) 



dU^ 
f7i 3) - TS„ 



-In Z< 3 > 



and 



S q = kln q 2f\ 

(3) 



In addition, use of Eq. fl27|), together with (3dU^> /d/3 = d(in q Z { q ] )/d(3 leads to 

5 



Note that Zg 3 ^ refers to the energy levels {ei}, with respect to Uq"* . We can choose as the energy reference 

(3) 

defining Z q ' through 
and recast (p7|) and (|29|) as, 

and 



(3) 



(26) 
(27) 

(28) 
(29) 



q ^q 

-in Z® 

m q ^q ' 



7j(3) 

u q 



dln g zj 3) 
9/3 



(30) 
(31) 

(32) 
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Finally, it can be shown that the Legendre structure is preserved too. The appealing features of this 
third choice are based in that it avoids all three problems that were noted for the second set of possible 
constraints pi]]. It is only necessary to define the new expectation values as 



Of ] - «Oi» 9 - (33) 



where O is any observable. Now, ((l)) q = 1 Vg and 

U® (A + B)= U® (A) + U® (B) . (34) 
It is important to mention that the probabilities associated with the third choice coincide with those 



obtained with the second if we use a normalized temperature [ 21 ] . This is the reason by which all theorems 



that do not use an explicit temperature dependence of a given phenomenon will continue to be valid. 
2.3 Distribution functions 

Despite that the third option is conceptually simpler than the second, actual computations of thermal 
dependencies are much harder. This stems from the fact that in the second option, the equations for {p^} 
are explicit (see fllED) whereas they are implicit in the third case (see (24)). The objective of this work is to 



analyze the evolution of the universe when the distribution functions slightly differ from the standard ones. 
To this end, working in the second option will provide the simplest expressions possible for all observables, 
what will allow analytical computations thorough. For particular cases, in addition, it has been proven 
that the difference between the exact results (using the third choice) and the approximation we shall make 
here are in very good agreement [^] . 

Let us now consider a gas composed by ./V non-interacting particles. If the system is in thermal contact 
with a heat and a particle reservoir, the energy and the number of particles will be conserved. The 
stationary state of the system will be given by the solution of 

HiP R = E R ^ R , (35) 

where H is the Hamiltonian and ipR is the wave function of the system. The accessible states will 
be represented by R. In second quantization formalism, R is given by the set of occupation numbers 
{ni, tt-2, . . . , Wjfc, . . .}, where denotes the number of particles in the state k. Taking n& = 0, 1 for fermions 
ensures Pauli's exclusion principle. We must now extremize the entropy with the imposed constraints, 

i = E p «' ( 36 ) 

R 

E = J2 P R E R> ( 37 ) 
R 

N = J2 P R N R- (38) 

R 

To obtain the solution we use the usual method, extremizing the expression (with units such that k = 1), 

-aY / PR-PT, P R E R-^T, P R N ^ ( 39 ) 
Ft Ft Ft 




where a,/3 are 7 Lagrange multipliers. Making, OQ/OPr = 0, we obtain 

E (^fj- + a + QP P R lp R + ^r'Nr) = 0. (40) 



R 
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Since it should be zero for all R, each term in the sum must vanish, 

qP^ 1 + (q - l)a + q(q - I^P^Er + 7 q(q - l)P Q R l N R = 0. (41) 
Identifying f3 and 7 by 

7 = -fin, (3 = 1/T, (42) 

where [x is the chemical potential and T the temperature, the probability that the ensamble is in the state 
R is, 

Pr = [1 + 0(q - 1)E R - P(q - l)fxN R f^ /Z q , (43) 

where 

Z q = £[1 + 0{q - l)E R - (3(q - 1) M ^]V(5-D. (44) 
R 

The occupation numbers automatically determine the quantum state of the system, 

Er = niei + n 2 e 2 + . . . + n k e k + (45) 
N R = m + n 2 + ... + n k + ..., (46) 

where e k is the energy of a particle in the state k. 

Using Eqs. (45) and (|46| ) in the expressions ( f44|) and (|43"|) we obtain, 

p ^ [1 + g(g - l)(ei - Am + ... + P(q-l)(e k - n)n k + . . 

with 



Z g = 2 [i + / g( g _i)( ei _ M ) ni + ... + / 3(g_i)( £fc _ M ) nfc + ...]i/(ff-i). (48) 

m,...,n kl ... 

The partition function can be factorized as 

00 

Z l=Il Etl + Z^-lXefc-zW^- (49) 

fc=l n fc =0 

It can be shown that (see Ref. [13] for details) the generalized occupation numbers are 

1 

with £ = 0, +1 or —1 in the case of a Maxwell-Boltzmann (MB), a Bose-Einstein(BE), or a Fermi-Dirac 
(FD) gas. In the limit, (q — 1) — ► standard distribution functions are recovered. The expressions for the 
distribution functions are not exact ones, since from Eq. ( |48| ) to (^) we have made and approximation, 
a.k.a. factorization approach |]i~3fl . In general, indeed, 

[1 + (q- 1)(A + B)] 1 ^-!) + [1 + ( 9 - l)^ 1 /^ 1 )^ + (g - l)^] 1 ^- 1 ). (51) 

Since the equality is valid for q = 1 we can certainly wait that for values \q — 1| <C 1 the approximation 
will remain a good one. We are interested in the first order in (q — 1), then the expressions (|50| ) becomes, 

~ e (3(e r -n) + £ + 2 (e/^r-M) + £)2 ^ 

In Fig. 1 we show the change in the distribution functions for different (q — l)-values. 
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-2 -1 1 2 -2 -1 1 2 0,0 0,5 1,0 1,5 2,0 2,5 

p(e-ji) p(e-n) P( E -H) 

Figure 1: Behavior of generalized distribution functions (FD, MB and BE respectively) for q = 1.1, 1.0 
and (q - 1) = 0.9. 



Eq. ( p2| ) is the expression we were looking for, it is simple enough, but yet accurate for small deviations, 
to give an interesting framework where to analyze (as a parameterization) the effect of a small change on 
the extensivity properties of the universe on its own evolution. Looking at Fig. 1 it can be seen that despite 
the particular analytical form for the new distribution functions, it is clear that they will do an excellent 
job in showing an slight deviation from the standard case, this being the ultimate reason why we use them. 

Up to recent days, quantum distribution functions within the factorization approximation were regarded 



as a rather rough technique, see Pennini et al. [23]. These latter authors considered fermion and boson 
systems with very small occupation numbers. However, Wang and Le Mehaute |^4| analyzed the problem 
in detail and showed that there exist a temperature interval, which they called 'forbidden zone', where the 
deviation from the exact result maybe significant. Otherwise, outside this zone, the factorization approach 
results can be used with confidence. In addition, they verified that the magnitude of the forbidden zone 
remained constant with the increase of the number of particles. This fact motivated new efforts in the 
study of macroscopic systems (where the number of particles is large), since the generalized distribution 



functions of the factorization approach could be used at temperatures up to 10 20 K for such systems [24] 
All these results encourage us further to use this approximation for the analytical study that follows. 



3 The plasma parameter, the plasma era, and before 

One of the earliest applications of a non-Maxwellian distribution to plasmas was referred to the solar core. 
Clayton considered an ionic quantum distribution function that slightly deviated from a BG one, being the 
correction factor proportional to exp[— 5(E/kT) 2 ] |^]. 5 can be related with a non-extensive parameter q 
by means of the formula S = (1 — q)/2, see Q. There are several reasons why to expect a slight deviation 
from a BG behavior in the sun core. Particularly, we know that long range forces are present (Coulomb 
and gravitational) and that there are also effects related with long range microscopic memory. Then, the 
leading idea is try to see whether the solar core is a non-ideal plasma, or a quasi-plasma, or equivalently, 
to what extent it can be considered as ideal. While we forward the reader looking for a detailed account 
of these effects to the references quoted in ||, let us just mention here the computation of the plasma 
parameter. 

The plasma parameter is the ratio between typical electrostatic energy of the nearest neighbors of a 
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given particle and the thermal energy kT. For a plasma with number density n, the mean inter-particle 
distance will be n 1 / 3 , and it can be proven that the number of particles inside a Debye sphere is given by 

m 



3/2 / „ X -1/2 



n °~*°[k) UfsJ • (53) 

The plasma parameter will be then defined as V = N d . When T <C 1, the plasma behaves (just 
considering the long range forces) as an ideal gas, the electrostatic energy is small compared with kT and 
there are a large number of particles that shield the long range force in a Debye-radius range. If, instead, 
r > 1, the plasma is strongly coupled and certainly non-ideal. For values of T near 1 we expect slight 
corrections. For the sun core, for instance, the plasma parameter is of the order of 0.1 Is it small enough 
to produce no appreciable deviation? The answer seems to be no. It have been recently proven, see last 
reference in Q, that the effects of random electric microfields are of crucial importance. These fields have 
in general long-time and long-range interactions that can generate anomalous diffusion that alters the 
distribution. Most importantly, the amount of the deviation has been quantified to be proportional to V 2 . 

The plasma parameter in the early universe is used to define what is known as the plasma era. At the 
end of the plasma era is where the process of recombination occurs. Here, the plasma parameter is very 
small and an ideal gas approach is, a priori, a very good one. If the universe is well described within the 
standard model of cosmology, we do not expect the distribution functions to deviate from the usual BG ones 
within the plasma era and in later epochs. Then, if we force the distributions to deviate, the corrections 
should turn out to be very big, in such a way that experiments could then impose restrictive bounds on, 
for instance, the non-extensive parameter. This is indeed what happens, as we shall see below particularly 
when analyzing the recombination process. Do the same happens for all eras of cosmic evolution? 

Near the temperature of freezing at about 1 MeV, something very important happens in the universe: 
the electron-positron annihilation. At about 0.5 MeV, the cosmological plasma changes its composition, 
before the annihilation being formed mainly by electrons, positrons, photons, neutrinos, and traces of non- 
relativistic particles. Which is the value of the plasma parameter in this era? This plasma is relativistic, 
and then we cannot expect that the previously introduced parameter (although indeed it grows as the 
temperature of the universe increases) to be a complete description. We can, however, consider it as a first 
approximation. Its value is ~ 0.07 at the time of freezing (t ~ Is, T ~1 MeV). This value of T is close to 
the sun core case. Can we a priori expect a similar correction to the quantum distribution functions? As 
we shall later show, when relativistic particles dominate the problem, corrections introduced by a slight 
deviation in the distribution functions are not so big as in the non-relativistic case. Restrictions on the 
parameter space are then not so strong, and, a posteriori, it is not possible to straightforwardly discard 
such slight deviations. In addition, if T 2 gives the amount of the deviation also for this system, we would 
expect corrections of the order of 10~ 2 — 10 -3 , what we do obtain (independently of these comments) when 
we analyze the epoch of primordial nucleosynthesis, see paper II in this series [ pTj] , 

We do not claim here that these paragraphs provide evidence of such an strong nature as to consider a 
priori that deviations to BG distributions functions are unavoidable, as the analysis for the sun core case 
seem to show ||. We do believe, however, that they make worth looking at the standard model predictions 
in an statistically modified setting, to see the allowed range of the deviations along the different cosmic 
eras. To that aim we devote this and the follow up paper in the series. On the other hand, we recall 
that several authors, with various motivations, have considered other effects that could slightly change the 
distribution functions as well: Non-equilibrium effects originating in residual interactions between neutrinos 
and electrons (see for instance [28]) and fluctuations in magnetic fields possibly present at that time (see 
for instance |29|) are some examples. 
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4 Thermodynamics in the expanding universe 

4.1 Energy density, pressure, number of particles 

The number density, the energy density, and the pressure for a dilute gas are given in terms of their 
corresponding distribution function /(p) by (units are such that h = c = 1), 



9 

n 



(27T) 



(2yr) 3 / fiP)d 3 P, (54) 
9 J f(p)E(\p\)d 3 p, (55) 



|2 



g 2 

(2*f c j 

where -E(|p|) = |p| 2 + m 2 , and g is the degeneracy factor of the particle. From now on, -E(|p|) = E(p) 
and |p| = p. Let us first take the non-degenerate case, in which kT 3> fj,. As before, T will stand for 
the temperature and fi for the chemical potential. Use of the form of f(p) in the non-extensive approach 
we adopted gives the correction to the standard result at first order in (q — 1). Several cases are worth 
mentioning. 

When the particles are relativistic, E » m, the energy density is 

1 5! 

Pq = Pst + ^2 -^( 1M 9bosons + 0.97g fermions ) (q - 1)T 4 , (57) 



with p st representing the usual result, 



Pst = ^ 9 T 4 (58) 



and g = gbosons + 7/8<7j erm j ons . For the number of particles, we obtain, 

n bosons = _^_ 2C(3)T 3 + ^12.98(g - 1)T 3 , (59) 

n fermions = _^_ 2 ^C(3)T 3 + _^H.36(g - 1)T 3 . (60) 

For the pressure, starting from the definition, we can immediately prove that P = p/3 disregarding the 
degree of non-extensivity. 

Remark 1: Note that, for relativistic particles, the temperature dependence of p,P and n in the non- 
extensive formalism is the same as the standard. This is because it comes only as a product of the change 
of variables in the integrals, which is done exactly in the same way for both contributions. 

It is interesting now to compute the energy per particle, which results to be 

< E q > boso ™ s = [2.70 + 11.29(g - 1)]T, (61) 

< E q > fermions = [3^5 + 4433^ _ !)] T . (62) 

Note that to equal q, the fermions get more affected than the bosons by the change of the distribution. 

In the non-relativistic limit m » T, we can neglect the ±1 term in the denominator of all integrals. 
For instance, the number density becomes, 
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9 



(mT\ 3/2 



e -(m-At)/T x 



1 + 



15 m — a 



m — p, 



(63) 



Here, the first term in the right hand side stands for the standard result and the second three for the 
correction at (q — 1) order. The energy density is then given by p q = mn q and the pressure is P q = n q T. 

Remark 2: For non-relativistic particles, the change to a non-extensive setting does change the tempera- 
ture dependence of the observables. This is not a minor effect and will show its impact in the recombination 
process. 



4.2 Particle anti-particle excess 

We focus now on the computation of the particle anti-particle excess. Consider the reaction r + +r~ «-> 7+7, 
where r + and r~ represent a generic particle and its anti-particle. If this reaction is in equilibrium, then 
the chemical potential of r + and r~ are related by p + = —pT . Then, the net number of fermions is given 
by 

2\ V2 . 



n' 



9 rp3 



2tt 2 



X \ x — — 



m/T 



T 



f+(x)-f-(x) 



dx, 



(64) 



where the integration variable is, as before, x = E/T. In the relativistic limit, m/T <C 1, and the standard 
integrals can be solved as 



■li- 



st 



n 



St 



6vr 2 



7T 2 — + ( — 

T \T 



To compute the (q — 1) order corrections we need to calculate 



9 T 3g~ 
2vr 2 2 



x 2 {x — if)) 2 e x ^ x 2 (x + tp) 2 e x+ ^ 



( e x--<p + iy 



dx. 



(65) 



(66) 



Since we would like to have an analytical solution, we think in a new application of the free parameter 
trick introduced in this context by Tirnakh and Torres j2^|. We know that 



j(m) 



2 m(«-i/>) _|_ \ 



3 m(u+t/)) _|_ 



(67) 



Here, m is considered as a free parameter. Then, the integral we need can be obtained by making the 
following steps. Consider —d 2 I m /dm 2 , evaluated at m = 1. This last integral will give the one we want 
plus two additional terms. These terms are 



du 



-2u 2 (u - ^) 2 e m («-V>) -2u 2 (u + ^) 2 e m ( u +^) 



( e m{u-i)) _j_ 1)3 



( e m{u+i>) _|_ 1^3 



(68) 



We can evaluate the relative strength of these extra terms with respect to the integrals we are looking 
for. The first thing one can analytically do is to consider the case in which ip ~ 0. This, we know, is 
not far from the real situation, since after electron positron annihilation, only a little excess of electrons 
will survive, in order to give the universe, together with protons, electrical neutrality. Because the ratio 
of the number density of protons to the number density of photons, n p st /n^ ~ 10 -9 , use of the standard 
expression for n 



st 



n st and n s J will show that 



n —^- ~ 1.33£ ~ 10- 9 , 
n^y T 



(69) 
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Figure 2: Relative strength of the integrals appearing in the computation of n+ — n c 



and clearly the condition fi/T <C 1 is sustained. Taking this into account we can weight each of these two 
extra terms with respect to our integrals. These terms are, 



and 



A 



B 



u A e u 



o (e u + l) 2 



Jo 



u A e u 



du, 



7 du. 



(e u + l) 3 

If we call /a(u) and /b(u) to the integrands of the previous expressions, we note that 

2 



f B {u) 



e u + 1 



fA(u). 



(70) 
(71) 

(72) 



Around 0, /s(u) ^ fA{u); but the integrand B falls exponentially. Finally, a numerical integration of the 
two terms, that we show in Fig. 1, yields a difference bigger than one order of magnitude: A = 22.63 and 
B = 1.037. For the sake of analytical computation, we shall adopt the following approximation for the 
correction to the particle anti-particle excess n c , 



c c ~ 2^ 2 



2 Tin 



-d Z I 



dm 2 



m=l 



which finally yields the result, 



^(1 + 3(1-,))^ + ^ 
In the case of non-relativistic particles, a similar computation gives 



n' 



( n tt ~ n st) + ( n t 



(73) 



(74) 



(75) 
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with the standard term being 



n st ~n st = g 



\ 2vr J 



sinh^/T), 



and the correction being 

'mT\ 3/2 



nt-n c =2g[—j 



-m/T 



sinh 



15 m / m? + u? 

T +3 T + hf^ 



cosh 



3^ + 2^ 
T T 



(76) 



• (77) 



4.3 Effective number of degrees of freedom 

The total energy density of all species in equilibrium can be expressed in terms of the photon temperature. 
As most of the energy density is contributed by relativistic species, we shall firstly consider only them. 
The total density will be the sum over all particular species, 



k bosons fermions 



(78) 



Note that we are using a sub- index i for boson particles, whereas j is used for fermions. From previous 
expressions we have, taking x r = E r /T, 



9rT, 



'4 ^oo 

/ X qfq {%r)dx r 

J m r IT r 



lTT z Jm r /T, 



(79) 



recall here that fg(x r ) is the non-extensive distribution function, and we are in the limit in which E r » m r 
and x r S> m r /T r are valid. If we as well consider T r » p r , then fg(x r ) = fg(x). This gives the separate 
corrections, 

Pbosons Pbosons, st + n _2 o" 

-5!C(5)7; 4 , 



2vr 2 2 
9 j Q — 1 

P fermions P fermions, st T" 2 

The standard contributions are the usual, 



5!/ 5 (l)T/. 



Pbosons, st 



30 



j ) P fermions, st — g - 1 j • 



Adding up both contributions, we can write 



(80) 
(81) 

(82) 



Pi? = ^9^T 4 , 



where gl is given by the addition of the standard, 



9* = 9i 

bosons 



Ti 



7 \- 

+ 8 ^ * 

fermions 



T 



and the new correction, 



5* = (9-1) 



9-58 £ +8.98 £ 

bosons fermions 



(83) 



(84) 



(85) 
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Remark 3: Note that, due to the same temperature dependence, we can define an effective number of 
degrees of freedom exactly in the same way as it is done in the standard case. Then, we can consider the 
same evolution equations, and hide all non-extensivity effects in the new gi. 



As an example of the previous remark, we can consider the Friedmann equation. It reads, 



R k 8vrG 
+ 



R R 2 3 



-P, 



(86) 



where R(t) is the scale factor in the FRW metric, and k is the curvature of the space-time. In the early 
universe, it is a very good approximation to take k = 0. During the radiation dominated phase P = p/3, 
R(t) = t 1 / 2 , and R/R = H(t) = 1/t. Using the expression for the energy density, i.e. Eq. (|83|), we obtain 



implying, 



H = l.m{glf/ 2 
t = 0.30 (g% 



mpi 

1/2 m Pl 
J>2 ' 



(87) 



Then, although the functional form is the same as in the standard case, the actual time-temperature 
relationship is different because of the change in the quantum distributions. 

5 Conserved quantities 

From the full set of Einstein field equations, or from the fact that T^ v . v = 0, it may be established that 

d( Pq R 3 ) = -P q d(R 3 ). 



This can be written as 



dT 



( Pq + P q )R 3 



R 



3 dPg 
dT 



And using the previous equation we finally get 



d_ 

dT 



Pq + P 1 R 3 



0. 



(89) 
(90) 

(91) 



as can be checked by direct differentiation. 



Remark 4-' Disregarding the degree of non-extensivity, s q = {p q + P q )/T is a conserved quantity in a 
comoving volume. 



If the relativistic contribution is dominant, 



bosons % fermions 3 



Pj 



This yields 

where we have defined gl )S as 



27^ 

"45~ 



q T 3 



9* s 9* s 9* s> 



(92) 

(93) 
(94) 
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with 



9*,st 



J2 9i 

bosons 



3 7 

+ 8 ^ * 

J ermions 



r n \ 3 

T 



g* = 7.18(g - 1) 



bosons 



TA 3 15 
+ 16 



71 \ 3 



(95) 



(96) 



/ermions 

From the conserved comoving quantity, we may define the number S q = s q R 3 = constant, and obtain the 
temperature-scale factor relationship as 

roc^j-^ir 1 . (97) 

Then, as in the standard case, T oc BT 1 , but now only in the periods where gl, s (and not <7* iS ) is constant. 
5.1 Explicit form for the conserved numbers 

For reasons that clarify themselves when actually making computations, it is convenient to define g% s such 
that g+ s = (q — l)p* s . As in the standard case, we have proven that for any particle which is not created 
nor destroyed during the evolution of the universe, 



Sq 



(98) 



is a conserved quantity. For future use, we quote here the explicit form of the conserved number for rela- 
tivistic particles. For instance, for photons, the number density is ri^ = g/{2-K 2 ) [2£(3) + 12£(4)(g — 1)] T 3 , 
with g = 2. Then, 

Sg = ^[2C(3) + 12C(4)( (7 -l)]- 1 ^X. 



This gives, 



N„ 



6C(4) - %(3) 

9*3 



with 

^ st ~ 2vr 4> 

Similar computations, but now for fermions, would yield, 



^ = ^-iC(% 6 (ri*)" 1 . 



with 



iv/« = ^W( 9 S)"- 



21CW - %C(3) 

9*s 



2vr 4 4 



(99) 
(100) 
(101) 

(102) 
(103) 



6 Decoupling 



Once a given species of particles has decoupled completely, each particle will travel along a geodesic of 
space-time. This will ensure that the form of the distribution function will remain the same during the 
subsequent evolution of the universe. Indeed, the only change will be given by a redshift correction. If 
the decoupling is produced at t = t d , then, due to the redshift, all particles that at the instant t have 
momentum p, should have had momentum pR(t)/R(td) at decoupling. Then, 

fdec(l>,t) = f eq {pR(t)/R(t d ),t d ), Vt>t d . (104) 

Remark 5: The previous result is independent of the particular f e q(p,t) we have adopted. 
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6.1 Relativistic decoupling 

Let us write the rhs of Eq. (|104| ) taken into account that E ~ p and in the case in which fi/T ~ 0, 

2 



p R(t)\ ; 



+ 



i 



(105) 



We see that the distribution function has the same functional form than f^ q but with a different "temper- 
ature" given by, 

R{t d ) 



T{t) = T d 



R(t)> 



(106) 



Despite this species is no longer in thermodynamical equilibrium, the "temperature" in the distribution 
function falls as -R -1 ; making Sg = s^R 3 to conserve separately. Note that for the species still in equilibrium, 
the temperature falls as T oc (<7* s (T)) _1 / 3 i? -1 . The number density for these particles is then given by, 



|- I [2C(3) + 12C(4)(g-l)]T 3 , 



IV, 



2vr 2 



~C(3) + §C(4)(g - i) 



rp6 



and using that for any time t, after t d , T(t) = T d , the previous equations really are, 



9 



^ 2 [2C(3) + 12C(4)( g -l)]T|(^ 



n 



fr-l 



fr 



2tt 2 



|c(3) + ^c(4)(«-i; 



T 3 / ^fa) 
" I i?(t) 



(107) 
(108) 

(109) 
(110) 



These number densities are comparable to that of photons at any given time. Such population will then 
continue to exist as a relativistic relic. 



6.2 Non-relativistic decoupling 

In the case of non-relativistic decoupling, the distribution function is given by 



fa(P,T)=e 



-{m-jj, d )/T d 



exp 



-P 2 ( R(t) 



2mT d \R(t 



1 + 



1 



p 2 ( R(t) 



2mT d \R{t a 



+ 



(m - fi d ) 



and the number density is given by 
3/2 /t>U w 3 



n. 



fmTdY 2 (R(t d ) 



(m-Hd)/T d 



1 + q ~ 1 (— + o m - ( m-fi d \' 

2 I 4 + T d + \ T d ) 



111) 



If we now consider that for a particle species which has already decoupled, the number density must 
be such that n q oc i? -3 , we see that as in the standard case (see for instance Section 3.4 of Ref. |30| ) , 
we need to identify the temperature with T d R 2 {t d )/ R 2 (t), and ask for the chemical potential to vary as 
fi = m — (m — fi d )T/T d . Then, the distribution function is similar to an equilibrium distribution but with 
a new temperature given by 



(113) 
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decreasing with the square of the scale factor. In the limit m 3> T, the energy density will be given by 
p = nm. 



Also in a non-extensive framework, the correct way to handle with decoupling process is through the 
Boltzmann equation, which we analyze below. 

6.3 The relic neutrino background 

The conservation of s q , applied to the particles which are in equilibrium with radiation shows that the 
quantity gi :S T 3 R 3 is constant throughout the expansion. During pair annihilation, g* jS will decrease, since 
we go towards a state with less number of degrees of freedom. Then, T 3 R? after the process will be bigger 
than its previous value. We say that photons are heated up by the annihilation process, in such a way that 

(gl s T*R 3 ) before = (gl s T^R 3 ) after (114) 

is sustained. This directly yields a relationship between both temperatures, 

rpafter fH^ 1 / 3 



T 



before 



-\ +0.152(g-l) = (-) (1 + 0.10%-!)) (115) 



Note that the first term is the standard result, from the usual hot big bang model. The second term is the 
non-extensive correction. 

Since before e + e~ annihilation, neutrinos were in thermal equilibrium with photons, T^ e f ore will be 
the same as j^ e / ore . Now notice that neutrinos do not participate in the process e + + e~ <-> 2j, and so, 
neutrino temperature is constant. Then, j'* e / ore = T^ ter , and we can write 



rpafter 



ll\l/3 

— +0.152(g-l) 



T^ ter . (116) 



After annihilation, g* ;S does not change anymore, and both temperatures T^ ter and T^ ter fall as R 1 , 
conserving their ratio. Since now T 7 ~ 2.728, then T u = 1.947(1 — 0.109((7 — 1))K. 

Remark 6: Even if the evolution of the universe diminish the degree of non-extensivity, the relic back- 
grounds, produced at early times, are sensitive to the possible non-extensivity present at that epoch. 

Remark 7: Note that if (q — 1) could have any value, one can have a photon temperature unchanged 
after e + e~ annihilation: T^ ter /j^ e / ore = i. This effect particularly show how important is the statistical 
description in the evolution of the early universe. Note, however, that a posteriori constraints on the value 
°f ( a ~ 1) do apply, and this situation will turn out to be impossible. 

The difference between T u and T 7 can be bounded using information coming from the CMBR, see for 



instance the work by Torres [11] 



7 Recombination 

We focus now on the generalization of the Saha law. Let n ,n p and n e the number density of hydrogen 
atoms, protons, and electrons respectively. Electrical neutrality implies that n e = n p . Baryon total number 
is given by n B = n H + n p . In thermal equilibrium, a temperatures such that T < rrii, we have seen that 
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is given by Eq. (p3|). The process p + e — > H + p guarantees that \i v -\- \i e = hh- We can as well consider 
that m p ~ rriH- Then, we would like to find an expression for n? jn e q n v . Immediate algebra yields, 



n 



H 



9H fm e 



n e q rf q g e g p \ 2ir 



exp 



m p + m e — ran 



T 



u((m H - Hh)/T) 



u((m p - n p )/T)u{{m e - u e )/T) 



Here, we have made use of the definition, 

Ui = u((m» - m)/T) = 1 + 



- 1 [15 „ Wi - jij 
2 I 4 T 



mi - m 



(117) 



(118) 



We would like to introduce now the number of baryons, since this is a conserved number (as far as precision 
measurements can tell), i.e. N B = n B /s is constant. Note however that r] q = n B /n^ does not remain 
unchanged, since gl tS may change in time. Using g p = g e = 2 and gn = 4 we can write, 



n 



nf, 



H 



B „B n 7^ 



n q 



2tt 



exp 



m p + m e 



u((m H - Hh)/T) 



u((m p - n p ) /T)u((m e - jJL e )/T) 



(119) 



We now define X q = nE/n B and 1 — X q = n B /n B Using them, we have 



li- 



st 



1 + 



q I '"q 

nL 



m e T\- 3 / 2 



B/T 



u n u, 



p <*e 



(120) 



where B = m p + m e — win and = n 7iSi + {q — l)n 7)C . Note that we have shortened the notation on the 
last bracket, but we are actually implying the definition of Eq. (|118j). Also, 



Vq = Vst 



1 + 



n 



B 



li- 



st 



1 + 



nj. 
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st 



Finally, we rewrite 



with 



and we use that 



1ip1l e 
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H T p~ 
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4 T ~ 



rrn - Hi 



(fiH—rn H )/T 
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g H eU*H-m H )/T _|_ g pe (p. P -m p )/T 



(121) 

(122) 
(123) 
(124) 



Recalling that nj t = (2£(3)/7r 2 )T 3 , we can write, with gn =, and g p = g e = 2, 



1 -xi 



{X e qY 



4^/2 



7T 



C(3)r?. 



'st 
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And taking into account the standard result, 
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.(125) 



(126) 
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Figure 3: Left: Standard fraction X^ t used to compute the corrections in the Saha law. It is worth noticing 
how quickly recombination is accomplished. This figure stands for a standard result, and it is quoted here 
just for the ease of the discussion. Right: Evolution of the standard chemical potentials as a function of 
the temperature for different values of Qh 2 . As far as we know, although this figure is an standard result, 
it was not presented in cosmology textbooks before. The numerical results contained in both of these plots 
were used to compute non-extensive corrections in the Saha law. 



we see that we have arrived to the following relationship: 

JxW + ~V^ Ci )Vst W 6 { 2e fr*-"*)/r + e (/»-™,)/r +VH-(u e + u p )).(127) 

This is the generalized Saha law. 

Looking at Eq. ( |127| ) we may note that we are going to confront a problem that does not appear in 
the standard case: The chemical potentials appear in a separate way, in such a way that not only the 
relationship hh — fj, e + fi p = enters in the computation. Since it is not common to find the ^-values, or 
even the method to find them, in standard cosmology books, we shall briefly outline here how they can 
be naively obtained. We shall then obtain the standard chemical potentials [/, = fi s t, and consider that 
generalized expressions for them will be written as [i q = fx s t + (q — 1)C. Since the term including the 
factor C will provide a second order correction in Eq. ( |127|) we shall be concerned only with the standard 
computation. We have to consider a system of three equations: 

1 . chemical equilibrium [in — \i e + fi p = 

2. electrical neutrality n e = n p 
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Table 1: Standard chemical potentials used in the determination of the correction at recombination epoch. 
Two particular values for Qh 2 are shown. 



T[eV] 


fi e —m e 
T 


fj,p Trip 
T 


T 


X? 


nh 2 


0.55 


-9.7450 


-21.018 


-6.0360 


1. 


0.1 


0.50 


-12.652 


-25.820 


-13.168 


1. 




0.45 


-20.434 


-31.707 


-21.918 


0.9998 




0.40 


-27.808 


-39.082 


-32.890 


0.9959 




0.37 


-33.156 


-44.429 


-40.829 


0.9481 




0.35 


-37.044 


-48.317 


-46.504 


0.7537 




0.33 


-40.991 


-44.429 


-40.829 


0.9481 




0.32 


-43.011 


-54.284 


-54.795 


0.2305 




0.31 


-45.107 


-56.380 


-57.616 


0.1267 




0.30 


-47.310 


-58.583 


-60.560 


0.0647 




0.25 


-60.806 


-72.077 


-78.486 


0.0008 




0.55 


-12.047 


-23.320 


-10.641 


1. 


i. 


0.50 


-16.850 


-28.123 


-17.773 


0.9999 




0.45 


-22.735 


-34.008 


-26.522 


0.9989 




0.40 


-30.076 


-41.349 


-37.426 


0.9619 




0.37 


-35.169 


-46.442 


-44.854 


0.7096 




0.35 


-38.658 


-49.931 


-49.732 


0.3786 




0.33 


-42.307 


-53.580 


-54.675 


0.1431 




0.32 


-44.251 


-55.525 


-57.277 


0.0797 




0.31 


-46.304 


-57.577 


-60.011 


0.0419 




0.30 


-48.484 


-59.757 


-62.908 


0.0209 




0.25 


-61.958 


-73.231 


-80.789 


0.0002 




0.20 


-82.191 


-93.464 


-107.655 


0. 





3. the standard result for 1 — Xg t /(X^ t ) 2 . The fraction X e st is shown in Fig. 3. 

These three equations form a complicated system but with three unknowns, regarded here as (/ij — rrn)/T, 
where i runs over the three species involved. Although this system cannot be solved analytically, it can be 
solved in a numerical way. The exact values of the solution will ultimately depend on the value of rj, the 
baryon to photon ratio, which in turns depends on the parameter VLh 2 . Results for the chemical potentials 
are shown in Table [l| and in Fig. 3. 



With these values, the correction terms in the Saha law, using Eq. (127), can be immediately computed. 
The results are shown in Fig. 4, where two different cases for the baryon content of the universe are 
studied. As the corrections in the non-relativistic expression for the number of particles is proportional to 
the factors ((/Xj — rrii)/T) 2 , we note that leading order corrections in the period we are interested in are 
~ (q — 1)50 2 = (q — 1)2500. These are large factors, and are the main reason by which all bounds imposed 
at the time of recombination, where matter and radiation cease to be in thermal contact, i.e. using CMBR 
data, are much stronger than those imposed at earlier epochs. 



8 The Boltzmann equation 

To analyze in an appropriate way the decoupling of particles we must consider the microscopic evolution of 
the distribution functions f{p^,x^) } via the Boltzmann equation. An interesting recent study on this issue 
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Figure 4: Final corrections in the Saha law. Note that the relative strength of the correction differs 
appreciably for different baryon densities. Note that the interesting range for the correction is limited by 
the equilibrium in the reaction D < — ► p + n, represented by fj,^ = /J-n + A*p- 

was presented in Ref. [ pL2[J , but their emphasis was put on other aspects. We shall start here by writing 
the Boltzmann equation as, 

L[f] = C[f], (128) 

where C is the collisional operator and L is the Liouville operator. The expression for the latter, in the 
non-relativistic case, for the phase space density f(v, x) of a species with mass m under a force F = dp/dt, 
is 

d dx -< dv -l , 

equivalently, 

d - F -. . . 

L = - + v-V s + --Vt. 130 
dt m 

The relativistic generalization becomes, 

Gravitational forces are present through the connection . For a Friedmann- Robertson- Walker universe, 
/ = f(\fi\it) ( or equivalently / = f(E,t)) and the only non-null components of the connection Fg are the 
usual 
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n> = f*,, (134) 



with h 13 = —g v and Latin indices running from 1 to 3. 



Taken this into account and recalling that p^ = (E,p) and = (t,x), Eq. ( 131 ) becomes 

L[f(E, t)\ = E^J- - -\p\ 2 ^J-. (135) 
uy ' n dt R m dE v ; 



Substituting this expression in Eq. ( 128 ) and taking into account the definition of the number density, 
we can integrate by parts to obtain 



" + :UIn = hl c[f] ^E- (136) 



We can also rewrite the collisional term defining 

f_J_C[f] , , 

in order to have the Boltzmann equation as, 

n + 3Hn = [ fd 3 p. (138) 



The collisional term is usually dominated by annihilations between particles and anti-particles {pp) and 
we assume now that the number of particles is identical to the number of anti-particles. In this way, f c is 
given by (see Section 9.2 of Ref. |3l| for details): 

fc = - I <°v> ffd 3 p, (139) 

where < av > is the product of the cross section times the velocity, averaged in velocities space. We can 
extract < av > out of the integral, evaluating it in an averaged energy. 

Let us focus in ^-particles and their anti-particles ■0. Taking into account Eq. (|139j) we integrate over 



momentum space in Eq. ( 138 ) to arrive at 

+ SHn^ = -<av> (nf) 2 . (140) 

Recall that the index q refers to the non-extensive distribution functions and generalized quantities. We 
can now account for the thermal production of particles. Naming T the term giving this contribution, we 
can write the Boltzmann equation as 

h\ + mn\ =-<av> (n^) 2 + T. (141) 

This term can be spelled out as in the standard case invoking thermodynamical equilibrium: in an 
static universe (R = H = 0), ri^ would be a constant and equal to n^ eq ; then, 

4 + 3Hnt = -<av> («) 2 - « e? ) 2 ) • (142) 

We can rewrite the previous equation in terms of the variable Y q = ri^/s q: where Y qm = n^ eq /s q , resulting, 

h% + Wr4 = -<av> S 2 q (Y 2 - Y q 2 eq ) . (143) 
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All this resembles the usual derivation, and indeed we can go still further. We can write + 
in terms of Y q . We shall take into account that s q R 3 = constant, and we derive with respect to time both 
the definition for Y q and the equation s q R 3 = constant. Finally recalling that H = R/R, we get 

+ 3Hn* = s q Y q . (144) 

Combining the last two expression for hf + 3Hn^, we obtain: 



dt 



<av> s q { Y„ 2 - Y} 



q.eq 



(145) 



Introducing explicitly the temperature, via x = m/T, with m being any appropriate mass scale, and 
using the radiation-dominated-era relation between x and T, Eq. (|88|), we get 

dt 



dx 



2(0.301)( 5 *) 



-1/2 "V , 



and we can write 



Y n 



dY q dx dY q 



q dt dt dx 
Using the Jacobian dt/dx, we can rephrase for Y q 



Y n 



H{m) dY q 
x dx 



where we have defined 



H{m) = 1.66(5 



q\X/2 



m 



rripi 



Note that H(m) is related with H by the simple expression 

H = X - 2 H(m). 

Finally, using Eqs. ( |145| ) and (|148| ) we can write the Boltzmann equation as 



dY, 
dx 



X < (TV > S q / 2 _ 

H{m) \ q q 



2 

q,eq) ) 



where < av > is given by 



< av >- 



V>\2 



fcd 3 p. 



(146) 
(147) 

(148) 
(149) 

(150) 

(151) 
(152) 



Considering other ipip annihilation channels, say tprp into a final state F, an additional term would appear 
that would be similar to that in the rhs of the previous equation, but with < cr^^ xS v > replaced by 
< cr^^pv >. Sum over all annihilation channels yields the final result in terms of the effective cross 
section < oav >, 



dY q _ x < ctav > s q 
dx H(m) 
Defining Ta = n %eq < a AV >, we get 



x dY q 

Yq e q dx 



H(x) 



•y2 -y^2 

1 1,eq 



Y 



Y 



eq , 



(153) 



(154) 



Remark 8: Note that this derivation is independent of the particular form for the distribution functions, 
but the final product has the same aspect than that obtained in the standard scheme, see for instance Section 
(5.2) ofRef. 
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8.1 Freezing 

Following the standard procedure, we would like to study how much the freezing times and temperatures 
are modified by a change in the statistical description. Using Eq. ( |154|) , we see that the comoving number 
density of tp particles is controlled by the ratio T/H, and the amount of the deviation of the distribution 
functions from their equilibrium values. When T/H is less than unity, the change in the number density 
is small, — AY q /Y q ~ — (xdY/dx)/Y qm ~ T/H < 1, annihilations stop and the number density freezes. 

The annihilation rate Ta varies as n q ^ eq times the averaged cross section < gav >• In the relativistic 
regime, n q ^ eq <x T 3 . In the non-relativistic regime, n q ^ eq oc e m < , and Ta exponentially decreases. In 
both regimes, Ta decreases with temperature, and eventually goes below the necessary rate to maintain 
equilibrium, which by definition occurs when x = Xf (the "freeze out"). Then we expect that for x < xj, 
Y q (x) ~ y gje g(x), whereas for x > Xf, Y q (x > x/) = Y q ^ eq {xf). 



We can compute the value of Y q ^ eq 
and Eqs. (|59|) and ([~~ 
We obtain: 



n 



q,eq 



/ s q in each of the previous cases recalling the definition of s q 



for relativistic {x <C 3) particles and Eq. (|63| ) for non-relativistic (x 3) particles. 

C(4). 



for relativistic bosons, 



q,eq 



Y f 

Q,eq 



0.2785° 



1 + 6 



C(3) 



(9-1) 



(g. 



0.278-gf(gl 



1 + 



,04) 
C(3) 



(9-1) 



(155) 



(156) 



for relativistic fermions, and 



1 q,eq 



0.1455x 3/2 e _x (5: 



1 + 



h 6x + x 

4 



(157) 



with x 3> 3, for non-relativistic particles. 

Since we can write g1 s as g1 s = g^ s + (q — l)g+ s , we can get a first order result for the above mentioned 
quantities, 



and 



bf„st\-i 



Y° = 0.278g°(g°l) 



Y[ eq = 0.278 i /(^*)- 1 



1 + (q - 1) 6 



C(4) 




C(3) 


git). 


R C(4) 




U C(3) 


9&) 



-ynr 
I q,eq 



V 2 p - X (n st )- 1 
*sJ 



O.U5gx s/2 e- x (g: 



1 + 



- 1 /15 9 3S, 

f — + 3x + x 2 - 

2 V 4 of* 



(158) 
(159) 

(160) 



with x 3> 3. 

In Fig. 5, we show the equilibrium abundance for non-relativistic particles together with the non- 
extensive correction. The correction diminishes exponentially with decreasing temperature, but this is an 
effect of the overall factor e~ x also present in the standard result, and not of the non-extensivity introduced. 



9 Current values 

We are interested in the value of the corrected ratio 



a, 



* with Pc = M, 
p c 87rG 



(161) 
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Figure 5: Standard equilibrium abundance and non-extensive correction for non-relativistic particles. 



where G = m~ z 2 with m pl = 1.2211 x 10 19 GeV and H(t) is H = 2.1332/i x 10" 42 GeV (0.4 < h < 1) 
From the current value of p^, we can get the contribution of relativistic particles 

8vrG 



3# 2 Hq 



1 q 



To get Pq we must know the the value of gi and g 1 

Carrying out a first order computation using Eqs. (|95"|), and (p6|), and Eqs. ( |S4"| ) and 
a universe populated by photons and three types of neutrinos, we get 



(162) 



i , and assuming 



(gDtoday = 3.36 [1 + 9.43(</ - 1)] (163) 

(9ls)today = 3.91 [1 + 7.11(g - 1)] (164) 

Taking the current value of the photon temperature, T = 2.728 ± 0.004 K, we can obtain: p^, s q , 
and Qq :R , 

(pf)today = 7.83 x 10- 34 [l + 9.43(g -1)] gem" 3 , (165) 

{s q )today = 2899.41 [l + 7.11(g-l)] cm" 3 , (166) 

{nDtoday = 411.77 [l + 5.4(g-l)] cm" 3 , (167) 

i^ q ,R)today = 4.17 x 10- 5 [1 + 9.43(g - 1)] /t" 2 . (168) 



10 Wimps in non-extensive statistics 



We now consider weakly interacting massive particles, for instance neutrinos with very small mass, in 
non-extensive statistics. 
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10.1 Relativistic relics, xj < 3 

In this case, Y qm is not changing with time, (see Eqs. ( |158|) and (|159| )), The asymptotic value of Y q , 
(Y q (x — > oo)) = y gi ooj is just the equilibrium value at freezing: 



Y q ,oo = Y q f eq (x f ) = 0.278-^ 



l + 8^(,-i; 



C(3) 



with 



x f < 3. 



The current abundance of the ip species today is 



n q,today ~ s q,todayY q ,oo, 



%,to day = 4^806.04 [1 + 14.31(9 - 1)] G&)" 1 . 



The mass density contribution of this kind of hot relics is given by 



Pq,today ~ n q,today m ip' 



and the fraction of the critical mass will be 



n 



q, today 



h 2 



8ttG 
3#n 



h 2 ^ 

" Vq^oday 



^q,todayh 



2 = jj/7.65 x 10- 2 [1 + 14.31(9 - 1)] (^) . 



(169) 

(170) 
(171) 

(172) 
(173) 



We know that QqIi 2 < 1; then we can apply this bound to the contribution to the species ip to Qoh 2 
and obtain a cosmological upper limit to the mass of this kind of particles 



< 13.08- [1 - 14.31(9 - 1)] gl s {g f T l . 



15, 

2 + -(2 + 2x3) 
lb 



Light neutrinos would decouple when T ~ MeV. At this temperature 

= 2 + 1 [2 + 2 + 2x3] + 

<&(*/) = 10.75 [l + 7.60( g -l)]. 
Then, for an species with two components, gf = g v = 2, and to first order in (q — 1), it results 

m v < 93.72 [1 - 6.7(9- l)]eV. 

This is a non-extensive generalization of the Cowsik and McClelland bound, see for instance Ref. p0| . 
Even when is improbable that an experiment with enough sensitivity could differentiate among values of 
9 using Eq. J177 ), further recalling that it is an upper bound what that equation is imposing, it is worth 
noticing that a change of statistics can have a direct influence upon the mass spectrum of the particles. 



(174) 

(175) 
(176) 

(177) 



10.2 Cold relics, x f > 3 

In this case, 



Y™Jx f ) = 0.145s* 3 / 2 e 



1 H — ( -r + 3x/ + xj 



(178) 



Now, Y qm strongly depends on m (since xj = m/T^). To make a numerical estimation we have to determine 
Tj, using the condition T ~ H. Reactions able to change the number of wimps A, are like AA <-> XX, 
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where X is a generic particle, assumed in thermal equilibrium. The average value of av can be expressed 
as 

< av >= a (—} . (179) 
\m J 

The value of the power law index k will depend on the details of the annihilation process, usually it is 
of order 1. The value of ctq depends on m and has a simple form for the extreme cases m <C mz and 
m S> mz, where mz — W 2 Gev is the mass of the Z-boson. 

For wimps with m < mz, the effective cross section o"o can be written as pq], 



(180) 



where Gf is the Fermi constant, G F = 292.8 GeV. The value of c depends on the fermion. We shall 
consider Dirac (spin 1/2) particles for which c ~ 5. 
The reaction rate is given by 

T = n q < av >, 

or, equivalently, using Eq. (|63|), 

/ m x 3/2-* 



p _ &09A T 3 
~ (2tt) 3 /2 



-m/T 



T J 



1 + 



1/15 m ( m\ 2 

~ t +3 t + T 



The expansion rate H as a function of temperature is 



(181) 
(182) 



# = i.66 (<?; 



Then, the condition T/H ~ 1 yields 



1 = 3.82 x 10- j 5a(<£) 



1/2-fc 



-m/T d 



1 + 



— 1/15 m 



a mm p i. 



(183) 



(184) 



Solving for e m / T[J , substituting the result in Eq. (|178| ) for and computing gq using Eq. (|180|) , we 

get 



(x f ) = 2.86 x lO- 9 ^)- 1 / 2 - 



q,eqi 

fe+1 / ™ \ -3 



/ m V 
VGevJ 

Given the current value s q ^today = 2899.41 x [1 + 7.11(q — 1)] cm" 3 , we obtain 



n„ 



— 6 «»- 1/2 (i) i+1 (^)^ 



1/15 m / m\' 
7.11 + - — + 3— + — 
2 I 4 ^ Td \Tj 



(185) 



(186) 



The energy of these particles today is p q> today = n q ^odaym, and the contribution to the critical density is 



^^q, today *l — o tj Pq,today 
one 







fe+i 



m \ ' 
GeV/ 



i + (g-i) 



1/15 m ( m\ 



(187) 



We would now like to solve for T d . Taking logarithm in Eq. (|184|), we obtain 



m 

— = 17.74 + 



In 



fJA 



+ 



1 



k In 



m 
T~ d 



+ 3 In 



m 



GeV 



+ In 



l + (<7-l) 



m 



1/15 m 
2\T +3 T d + \T d J 



(188) 
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Since gi is a slowly varying function we can solve this equation in an iterative form. Consider the case 
of a wimp with m > 1 GeV. From Eq. ( |l88|) we see that m/Td ~ 17.74, assuming for simplicity that 
k = 0, the term ln(m/Td) corrects the value m/Td to m/Td ~ 19.18, giving Tj ~ 52 MeV (m/GeV). At 
this temperature, all species are relativistic and all of them contribute to gi, making it of the order of 100. 
Taking this into account, the term ln((^) 1 / 2 Jqa) — In 5 ~ 1.61, and the value of m/Td gets corrected to 
m/Td ~ 17.57. The non-extensive term is 



hi 



! + (<?-!) 



1/15 m 
2 T + 3 XJ + 



Substituting this expression in the zero order of m/Td we get 



rn 



= 17.57 + 3 In 



m 
GeV 



+ (q- 1)185.8. 



For masses m ~ 1 GeV, 



m 



— = 17.57+ (q- 1)185.8. 

Td 



The use of this result in Eqs. ([185]) yields 

4.87 x 10~ 9 [1 + W.5(q - 1)] 



-ynr 
1 q,eq 



m 
GeV 



n 



q,today 



}? = 1.34 [1-75.2(5-1)] (^_)' 



(189) 

(190) 
(191) 

(192) 
(193) 



The fermion A and its anti-particle A, will provide twice this value of Qqjodayh 2 , i.e. the constraint 
Qq,todayh 2 < 1 is then a corrected lower bound for the mass of these particles 



m > 1.64 [1 - 37.6(<j - 1)] GeV. 



(194) 



11 Modification of the matter-radiation equality 

Let us characterize the matter-radiation equality by the time of its occurrence t = t eq , its scale-factor 
R = R eq , and its redshift z = z eq . From the equalities 



3 

^3 



Pmatter{teq) — Pmatter(to) \ f ) — Pc^matter (to) (1 + Z e<? ) a , (195) 

\K(t eq ) J 

(R(t ) \ 4 
— Pc 5 ^matter 
(t )(l + z eq )\ (196) 
^(leg) / 

it follows 

_ ^matter (*o) _ » , 

J L radiation\ lj § ) J L radiation\ L ) 

And then, 

1 + z eg = (n fl (to)/i 2 )-W. (198) 

In the non-extensive setting, 

1 + z eq = 2.4 x 10 4 [1 - 9.43(<? - 1)] Vth 2 (199) 
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corresponds to a temperature given by T eq = Tq(1 + z eq ). Taking To = 2.728 K, we obtain 

T eq = 5.64 [1 - 9.43(g - 1)] Qh 2 , 



(200) 



and for the time we get 



t eg ~ QMH^n' 1 / 2 ^ + z eq )~ 3/2 = 1.022 x 10 3 [1 + U.U{q - 1)] (ft/i 2 )" 2 years. 



(201) 



12 Conclusions 

In this paper we have studied how the statistical description that we suppose is valid in the early moments 
of the universe affects different processes along the thermal history of our basic cosmological model. In 
particular, starting from the modifications on the energy and number densities, we explored all consequences 
introduced by non-extensivity in the processes of decoupling, particle-anti-particle excess, recombination, 
photon-neutrino temperature relationship, matter-radiation equality, and others. 

In order to avoid re-stating our obtained results, let us just mention a general conclusion. Suppose that 
all cosmological observables are measured with a precision of 1%. Suppose, too, that the standard values 
for these observables agree well with observations, i.e. the standard values are always within the error bar of 
the experiments. (This latter assumption is not always true, particularly when more than one experiment is 
involved at the same time, like in the case of the determination of the baryon density using nucleosynthesis 
and recent CMBR measurements.) When relativistic particles are involved, we can summarize our results 
as follows. Given any standard prediction, the correction introduced by non-extensivity is of the order of 
ten times the standard value, 



where the factor F is order 10X s t. Then, given the previous assumptions, when relativistic particles are 
involved, \q — 1| is bound to be of the order 10 -3 . On the other hand, when non-relativistic particles are 
involved, the corrections are much bigger, and F can be at least of the order of 1000X, t . In these cases, 
and also within the previous assumptions, \q — 1| is bound to be less than 10~ 4 . It is worth noticing that 
this simple conclusion is sustained by all currently existing bounds, and it is also worth noticing that the 
ultimate reasons by which the constraints on \q — 1| are more restrictive in the epoch of recombination 
than those obtained in nucleosynthesis (see Paper II) can be tracked from this formal point of view. 

In our second paper in this series we shall study the nucleosynthesis epoch, and will provide an analyt- 
ical assessment of the Helium 4 and Deuterium primordial production in the framework of non-extensive 
statistics, including corrections coming from the free neutron decay process. We shall particularly focus 
on how non-extensivity affects the principle of detailed balance. 

Appendix: Useful mathematical formulae 

It is convenient to introduce the following functions: 



X = X st + (q-l)x F, 



(202) 





It is immediate to show that lim^i e q = e x and that lim^i ln ? x = lnx. Also the following properties arc 
valid 
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With these definitions, 



S q = kln q W. 



(206) 
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